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Abstract 

We present an ab-initio density function theory to investigate the electronic and magnetic struc- 
tures of the bilayer graphene with intercalated atoms C, N, and O. The intercalated atom although 
initially positioned at the middle site of the bilayer interval will finally be adsorbed to one graphene 
layer. Both N and O atoms favor the bridge site (i.e. above the carbon-carbon bonding of the lower 
graphene layer), while the C atom prefers the hollow site (i.e. just above a carbon atom of the 
lower graphene layer and simultaneously below the center of a carbon hexagon of the upper layer). 
Concerning the magnetic property, both C and N adatoms can induce itinerant Stoner magnetism 
by introducing extended or quasilocalized states around the Fermi level. Full spin polarization can 
be obtained in N-intercalated system and the magnetic moment mainly focuses on the N atom. In 
C-intercalated system, both the foreign C atom and some carbon atoms of the bilayer graphene are 
induced to be spin-polarized. N and O atoms can easily get electrons from carbon atoms of bilayer 
graphene, which leads to Fermi level shifting downward to valence band and thus producing the 
metallic behavior in bilayer graphene. 

PACS Numbers: 71.15.Mb, 73.43.Cd, 81.05. Uw 



Typeset by REVTp^X 



1 



Introduction: Since the discovery of monolayer graphene in the year 2004 [2J, this 
two-dimensional material remains in the focus of active research motivated by its novel 
physical properties and a promising potential for applications^ 5]. Experimental groups 
have enabled preparation and study of systems with one or a small number of graphene 
layers [6|. Bilayer graphene, which is made of two stacked graphene layers, is considered to 
)e particularly importance for electronics applications because of its special band structure 
3]. Coupling of the two monolayer graphene sheets in the usual A-B stacking of bilayer 
graphene yields pairs of hyperbolically dispersing valence and conduction bands that are 
split from one another by the interlayer interaction [s]. The band gap of the bilayer can 
be easily opened when a difference between the electrostatic potential of the two layer is 
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introduced, either by chemical doping or by applying gate voltage I9H17H. which makes 
graphene channels have a high resistance for the OFF state. Angle-resolved photoemission 
(ARPES) measurements indicate such a gap in potassium doped bilayer graphene epitaxially 
grown on SiC [jjj], and infrared spectroscopy measurements also detect the similar gap in 
the electrostatically gated bilayer graphene [10h12|. 

In this work, we carried out first-principle calculations and theoretical analysis to explore 
the electronic and magnetic properties of bilayer graphene with intercalated atoms C, N and 
O. The calculations were performed using the projector augmented wave (PAW) formalism 
of density functional theory (DFT) as implemented in the Vienna ab initio simulation pack- 
age (VASP) [ISf. Because generalized gradient approximation (GGA) [19[ gives essentially 
no bonding between graphene planes and leads to excessively large values of bilayer distance 

131 ] . we performed the calculations within the localized density approximation (LDA). We 
find that LDA gives rise to a bilayer distance of 3.34 A, in good agreement with the ex- 
periment value. Some previous investigations used LDA to optimize the structure to get a 
reasonable bilayer distance, and sequently used GGA to calculate the electronic structure 

9[. To keep consistent, we use LDA in all our calculations for the investigated system. An 
energy cutoff of 400 eV for plane-wave expansion of the PAWs is used. The model system 
here consists of a 4 x 4 supercell with a foreign atom intercalated between the two coupling 
graphene layers. The supercell parameters are set to be the same as a = b = 9.84 A in 
the xy plane (a and b indicate the crystal lattice constants). The Brillouin zone is sampled 
using a 11 x 11 x 1 T centered k-point grid. For geometry optimization, all the internal 
coordinates are relaxed until the Hellmann-Feynman forces are less than 0.01 eV/A, The 
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vacuum thickness along the z axis is 16 A to avoid the interaction between graphene layers 
of adjacent super cells. 



We considered the A-B Bernal stacking structure for bilayer graphene. Top-view for the 
bilayer graphene is shown in Fig. 1(a), in which violet color is for the lower layer, and gray 
color is for the upper layer. Seen from the top-view, 'A' position in the lower layer coincides 
with 'B' position in the upper layer, 'A' position in the upper layer is exactly the center of 
the hexagon of the lower graphene layer from the top-view, and 'B' position in the lower 
layer is right the center of the hexagon of the upper graphene layer. The foreign atom is 
intercalated between the upper and lower graphene layers and is initially positioned in the 
middle position of the bilayer distance. Three kinds of initial positions are considered for 
each intercalated atom: bridge, top and hollow positions, which are noted in the following by 
using subscript "Bri" , "Top" , and "Hoi" when needed, for example, we use N Br i to indicates 
the bilayer system with N atom positioned at the bridge site. Bridge site is above the 
carbon-carbon bonding in the lower graphene layer, top site indicates the middle position 
between the coinciding positions 'A' and 'B', and hollow site is just above a carbon atom 
of the lower graphene layer and simultaneously below the center of a carbon hexagon of the 
upper layer. The structure relaxation calculations show that, if the intercalated atom is 
initially positioned at the bridge or hollow site, it will finally be adsorbed to one graphene 
layer and away from the other layer. While for the top site, the foreign atom will keep being 
located in the middle position of the bilayer interval. 

The binding energy is defined as: dE = E graphene + E atom - E total , where £ grap h e ne is the 
energy of the clean bilayer graphene, E atom stands for the energy of the single foreign atom, 
and Etota\ is the total energy of the bilayer graphene with the intercalated atom. The binding 
energy, the distance between the foreign atom and its nearest C atom, and the interlayer 
distance of the doped bilayer graphene are illustrated in Table I. Comparing the data in 
Table I, we find both N and O atoms favor the bridge site, while the C atom prefers the 
hollow site. No intercalated atom is stable at the top site. At the top position, the relaxation 
calculations show that nearest C-N bonding is about 1.83 A, and C-0 bonding is about 1.74 
A, which is much larger than the typical lengths of C-N (1.47 v4)and C-0 (1.42 A) 20], 
implying the physisorption rather than the chemisorption. For the hollow site, O atom is 
adsorbed to the layer which provides the nearest C atom, and away from the hexagonal 
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cental of the other graphene layer. atom and the nearest C atom form the C=0 bonding 
with the length of 1.39 A, smaller than the the length of C-0 bonding in the Oen structure 
(1.44 A). The calculation shows N Ho i structure doesn't exist because of the negative binding 
energy. Concerning the bilayer distance, with the foreign atom intercalated in the bilayer 
space, the distance is enlarged for all the investigated systems. 

Figure 1 displays the ground states for C, N and O configurations, where Fig. 1(b), 
(c) and (d) are the Ch i, Neri, and Oeri structures, respectively. In the Ch i system, the 
optimized interlayer distance is about 3.39 A, which is a little larger than the value of 3.34 
A for the pure bilayer graphene. Calculations show that C-C bonding in the upper layer 
nearly keeps unchanged. The carbon atom which is right below the foreign C is pushed 
down, forming a "dumbbell" at the saddle point. The foreign C bonds with the adjacent 
three carbon atoms with bonding length being 1.54 A, which is a standard length for sp 3 
hybridization [20J. In Neri system, the bilayer distance is enlarged to the value of 3.87 A, 
and the length of C-N is 1.43 A, indicating the chemically adsorption not the physically 
adsorption. The two carbon atoms bonded with N atom are drawn out of the graphene 
layer, and the C-C bonding is 1.55 A, which implies sp 3 hybridization of the two carbon 
atoms. In Oe r i system, bilayer distance is about 3.76 A, length of C-0 bonding is 1.44 A, 
and the two carbon atoms bonded with O atom are also drawn out of the graphene layer 
with the C-C bonding is 1.51 A. Obviously, N and O atoms have the similar ground state 
structure. 

Side-view of the charge distributions for Ceri, Neri, and Osri systems are displayed in Fig. 
2. Although Ceri structure is not the ground state for C-intercalated system, we show its 
charge contour together with those of N Br i and Obti systems, to provide a clear comparison. 
We know O atom is lack of two electrons, when it is adsorbed in the bilayer graphene, it 
strongly interacts with its adjacent carbon atoms and get electrons from them (see the upper 
panel of Fig. 2). N atom lacks three electrons, it also gets electrons from its adjacent carbon 
atoms (see the middle panel of Fig.2). C atom lacks four electrons to get saturated, and it 
tends to share electrons with its adjacent carbon atoms in bilayer graphene (see the lower 
panel of Fig. 2). All these three atoms show strong interaction with lower graphene layer. 
Concerning their interactions with the upper graphene layer, judging from Fig. 2, C atom 
has the strongest interaction with the upper layer, N is the second, and O is the third. We 
compare the stable positions for C, N and O atoms in monolayer and bilayer graphene, as 
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illustrated in table II. It is clear that N atom has the same stable position in both monolayer 
and bilayer graphene, so does the O atom. C adatom is stable at the bridge site in monolayer 
graphene, and Hollow site in bilayer graphene. We believe the different stable positions for 
C atom results from its interaction with the upper graphene layer. 

The Choi structure is favored as the ground state for C atom intercalated in the bilayer 
graphene. The spin-resolved band structures and density of states (DOS) for Choi system 
are shown in Fig. 3. Fig. 3(a) and (b) are the band structures for the majority spin and 
minority spin, respectively. It is clearly seen that impurity bands for the majority and 
minority spin components lie, respectively, lower and higher than the Fermi level. In the 
spin- resolved total DOS (see Fig. 3(c)), two narrow peaks at the opposite side of the Fermi 
level are observed, indicating the itinerant magnetism triggered by the intercalated C atom. 
In the following, we will find the itinerant magnetism comes from not only the foreign C 
atom, but also the carbon atoms of the bilayer graphene. In the orbital-resolved PDOS 
of the foreign C atom, s state and p state have peaks in the same energy range, which is 
indicative of the sp 3 hybridization. The C-C bonding between the foreign C and the nearest 
carbon atoms with the value of 1.54 A, is also a proof of the sp 3 hybridization, which has 
been pointed in the previous analysis. When the spin degree of freedom is neglected, our 
calculation from first principle predicts a twofold degenerate peak at the Fermi level (Fig. 
3(e)), but the spin unpolarized state is not the ground state. Including the spin degree of 
freedom, the balance between the majority and minority spin components will be destroyed. 
The spin density distribution of the lower graphene layer of the Ch i system are shown in 
Fig. 3(f). Both the foreign C atom and the pushed down carbon atom are magnetic, yet 
we cannot see their magnetic moment distribution in Fig. 3(f), because they are not in the 
lower graphene plane. Seen from the top-view, these two carbon atoms occupy coinciding 
positions, which are noted in Fig. 3(f) by the letter 'C The three nearest carbon atoms 
bonded with the foreign C atom are nearly nonmagnetic and the next-near-neighbors are 
magnetic. On the whole, the total magnetic moment of the Choi system is about 1.32 fiB, 
and the magnetic moment distributions show threefold symmetry, which is similar with the 
hydrogen adsorption on the graphene plane j^J. 

Spin-resolved band structure and density of states for Neri system are shown in Fig. 4, 
where Fig. 4(a) and (b) show the band structure for majority and minority spin, respec- 
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tively, and Fig. 4(c) and (d) display the total DOS of Ne r i system and PDOS of N atom, 
respectively. In the total DOS, very narrow and sharp peak at the Fermi level is observed. 
We draw the PDOS of N atom and find the peak arises from from the N adsorption and 
the N atom is nearly 100% spin polarized. In the band structure, we find the very localized 
states at the Fermi level in the minority spin band structure. Such quasilocalized states give 
rise to the strong Stoner magnetism with magnetic moment of 0.65 /iB located at the N 
atom. In the PDOS of N atom, between the energy -0.5 eV and -1 eV, another narrow peak 
is obtained. Both in the majority and minority spin band structures, the corresponding flat 
impurity bands in the energy range from -0.5 eV to -1 eV are observed. In addition, both in 
the majority and minority band structures, the characteristic conical point at K point still 
can be clearly identified, implying the bilayer graphene is not strongly purturbed by the N 
atom. The critical difference is that, in freestanding graphene the Fermi level coincides with 
the conical point, while the Fermi level obtained in N-intercalated system is shifted down- 
ward and becomes below the conic point. A shift downwards (upwards) means the holes 
(electrons) are donated by the adsorption atom. The manganese doping (electron donor) 
results in the upward shift of the Fermi level, reported by previous investigations {9]. 

Band structure and density of states of Oeri system are displayed in Fig. 5, where Fig. 
5(a), (b) and (c) are the band structure, total density of states, and partial density of states 
of O atom, respectively. The calculation shows the Obo system is nonmagnetic. From the 
above analysis about Neri and Ch i systems, we see the impurity states localized around the 
Fermi level play an important role in the magnetic properties, and the intercalated atom 
itself is spin-polarized. In Oeri system, we do not get such localized states around the Fermi 
level, and the O atom is nonmagnetic. In the PDOS of O atom, we find the peak nearest 
to the Fermi level is around the energy of 0.75 eV. From the band structure, we can also 
see that the impurity state nearest to the Fermi level is at the energy of 0.75 eV. Like Neri 
system, the Obi-i system also undergoes a change from semimetal to metal, because O atom 
get electrons form graphene, thus the Fermi level shifts downward. 

In summary, we calculated the structure, electronic and magnetic properties of the bilayer 
graphene with intercalated atoms C, N, and O by the ab-initio density function theory. 
Impurities even at very low density may bring fruitful physical properties for bilayer graphene 
system. Structure relaxation shows that the intercalated atom although initially positioned 
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at the middle position of the bilayer distance, will finally be adsorbed to one graphene layer 
and away from the other layer. The bilayer distance is enlarged for all the investigated 
adsorption systems. Both N and O atoms are stable at the bridge site, while C atom 
is stable at the hollow site. Concerning the magnetic property, O-intercalated system is 
nonmagnetic, while N and C atoms can induce stoner magnetism by introducing extended 
or quasilocalized states around the Fermi level. Nearly 100% spin polarization is obtained in 
N-intercalated system and the magnetic moment focuses on the N atom. For C-intercalated 
system, the magnetic moment distributes on the foreign C atom and certain carbon atoms 
of the bilayer graphene. Both in N-intercalated and O-intercalated systems, the Fermi level 
is obviously shift downward, inducing the metallic behavior of the bilayer graphene. 
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I. FIGURE CAPTIONS 



Fig. 1 (a) A-B Bernal stacking structure for bilayer graphene. (b) Ch i system, where the 
atom in violet is the foreign C atom, (c) N Br i system, where the atom in blue is the 
foreign N atom, (d) Br i system, where the atom in red is the foreign O atom. 

Fig. 2 Side-view of charge contours for C Br i, N Bri and Bri systems, respectively. 

Fig. 3 Electronic and magnetic structures of the Ch i system, (a) Band structure of ma- 
jority spin, (b) band structure of minority spin, (c) total density of states, (d)the 
orbital-resolved DOS for the foreign C atom, (e) density of states for spin unpolarized 
Choi system, (f) distribution of the spin density on the lower graphene layer.. 

Fig. 4 Spin-resolved band structure and density of states of N Bri system, (a) Band structure 
of the majority spin; (b) band structure of the minority spin, (c) total density of states, 
(d) partial density of states of N atom. 

Fig. 5 Band structure and density of states of Bri system, (a) Band structure, (b) total 
density of states, (c) partial density of states of O atom. 

Table I The binding energy (dE), the length of bonding between the adsorption atom and 
its nearest C atom (ac-atom), and the bilayer distance of the adsorption system (dis.). 
The stable position for each intercalated atom is noted by box. 

Table II The stable positions for adsorption atoms C, N, and O in monolayer and bilayer 
graphene. 
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